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ABSTRACT 



Context. In some recent works, the C/0 abundance ratio in high-metallicity stars with planets is found to vary by more than a factor 
of two, i.e. from ~ 0.4 to C/0 ~ 1. This has led to discussions about the existence of terrestrial planets with a carbon-dominated 
composition that is very different from the composition of the Earth. 

Aims. The mentioned C/0 values were obtained by determining carbon abundances from high-excitation C i lines and oxygen abun- 
dances from the forbidden [O i] line at 6300 A. This weak line is, however, strongly affected by a nickel blend at high metallicities. 
Aiming for more precise C/0 ratios, oxygen abundances in this paper are derived from the high-excitation O i triplet at 7774 A. 
Methods. The Ci lines at 5052 and 5380 A in HARPS spectra were applied to determine carbon abundances of 33 solar-type stars 
for which FEROS spectra are available for determining oxygen abundances from the O i X1114 triplet. Differential abundances with 
respect to the Sun were derived from equivalent widths using MARCS model atmospheres. Non-LTE corrections were included, and 
the analysis was carried out with both spectroscopic and photometric estimates of stellar effective temperatures and surface gravities. 
Results. The results do not confirm the high C/O ratios previously found. C/0 shows a tight, slightly increasing dependence on metal- 
licity, i.e. from C/O 0.58 at [Fe/H] = 0.0 to C/O 0.70 at [Fe/H] = 0.4 with an rms scatter of only 0.06. 

Conclusions. Recent findings of C/0 ratios higher than 0.8 in high-metallicity stars seem to be spurious due to statistical errors in es- 
timating the strength of the weak [O i] line in the Ni i blended /t6300 feature. Assuming that the composition of a proto-planetary disk 
is the same as that of the host star, the C/0 values found in this paper lend no support to the existence of carbon-rich planets. The small 
scatter of C/O among thin-disk stars suggests that the nucleosynthesis products of Type II supemovae and low- to intermediate-mass 
stars are well mixed in the interstellar medium. 
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1. Introduction 

Recent high-resolution spectroscopic surveys suggest that the 
carbon-to-oxygen abundance ratio, C/CQ, in the atmospheres of 
metal-rich F and G dwarf stars varies by more than a factor of 
two. Delgado Mena et al. (120101 ) find a range in C/O from ~ 0.4 
to i 1 for 370 stars observed with HARPS or CORALIE, in- 
cluding about 100 stars with detected planets. A similar range in 
C/0 has been found by Petigura & Marcy (1201 11 1 for about 700 
stars observed with HIRES of which 80 have detected planets. 
Even if one corrects the distiibution of C/O for the quoted ob- 
servational errors, there seems to be a small fraction (1-5%) of 
the stars having C/O i 1, and a significant fraction (10 - 15 %) 
having C/O > 0.8 (Fortnev l2012l l in contrast to the solar ratio of 
C/0 ^ 0.55. 

Assuming that the composition of a proto-planetary disk is 
the same as that of the host star, variations in the C/O abundance 
ratio among solar-type stars may have interesting effects on the 
composition and structure of planets. Thus, Kuchner & Seager 
(120051 1 discuss the possible existence of extrasolar "carbon plan- 



* Based on data products from observations made with ESO 
Telescopes at the La Silla Paranal Observatory under programmes given 
in Tabled 

** Tables [T] and |4] are provided as online material and is available in 
electronic fom at http : / /www . aanda . org 

' C/0 is defined as Nc/Nq, where A'c and A^o are the number densities 
of carbon and oxygen nuclei, respectively. It should not be confused 
with the solar-normalized logarithmic ratio, [C/O] = log(A'c/A'o)stai - 
log(A'c/A'o)sun- 



ets" consisting of carbides and graphite instead of "Earth-like" 
silicates. When modelling terrestrial planet formation. Bond et 
al. (20 10) find that C/O > 0.8 leads to the formation of such 
carbon planets. A particularly interesting case is super-earth 55 
Cancri e. Studies of its transit across the host star (Winn et al. 
201 1; Demory et al.| 201 l| l lead to a radius R ^ 2.0/?® and amass 
M ^ 8.4 Me. Adopting an oxygen-rich composition, a mas- 
sive water envelope over an Earth-like inteiior is required to ex- 
plain the observed radius. Madhusudhan et al. (12012) have, how- 
ever, shown that the mass and radius can also be explained by 
a carbon-rich interior without a volatile envelope. Interestingly, 
the star hosting this planet, i.e. 55 Cancri alias HD75732, has 
C/O =1.12 according to Delgado Mena et al. ([2010D . 

The existence of dwaif stars with C/O > 1 has been ques- 
tioned by Fortney (I2012I I. Among cool K and M dwarfs, C/O > 
1 leads to strong C2 and CN bands that are easy to detect. 
From various large spectroscopic surveys, including the SDSS, 
Fortney notes that the relative frequency of such "carbon" stars 
is less than ~ lO""*, and discusses various reasons that the C/O 
ratio may have been overestimated in solar- type stars. 

In the works of Delgado Mena et al. ( 120101 1 and Petigura 
& Marcy (1201 11 1, carbon abundances were determined from 
high-excitation C i lines and oxygen abundances from the zero- 
excitation, forbidden [O i] line at 6300 A. The large difference in 
excitation potential makes the derived C/O sensitive to errors in 
stellar effective temperature and to the temperature structure of 
the model atmospheres. Although the analysis in the two works 
was made differentially to the Sun, the use of homogeneous ID 
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model atmospheres might have introduced systematic eiTors of 
C/O as a function of metallicity. Furthermore, the [O i] line is 
overlapped by a Ni i line, which is calculated to have a strength 
of about 55% of the [Oi] line in the solar flux spectrum (Caffau 
et al. l2008l l. if one adopts the gf/-value determined by Johansson 
et al. ( 120031) . Given that [O/Fe] decreases and [Ni/Fe] increases 
with increasing [Fe/H], (e.g. Bensby et al. 2005), the Nil blend 
becomes more important at the high metallicities for which most 
of the high C/O values have been found. This makes it difficult 
to determine precise oxygen abundances from the [O i] line at 
high metallicities, and systematic errors in [O/H] are introduced 
if the Ni i blend is not correctly calibrated for the solar spectrum. 
Finally, the measured equivalent width of the weak [O i] line de- 
pends critically on the continuum setting, and sometimes it is 
difficult to correct for overlapping telluric O2 lines. 

Instead of the [O i] line, the O i triplet at 7774 A may be used 
to determine oxygen abundances. These three lines have a high 
excitation potential like the C i lines, which makes the derived 
C/O ratio insensitive to errors in effective temperature and the 
temperature structure of the model atmospheres. Furthermore, 
the lines are practically unblended. A drawback is that the triplet 
is affected by deviations from local thermodynamic equilibrium 
(non-LTE), (Kiselman 2001), which is not the case for the for- 
bidden [Oi] line. 

The ESO 3.6m HA RPS sp ectra (Mayor et al. [20031 used 
by Delgado Mena et al. (120101 ) have a maximum wavelength of 
6900 A, but for a good fraction of their stars, ESO 2.2m FEROS 
spectra (Kaufer et al. 1999) covering the Oi triplet are avail- 
able. In this paper, these FEROS spectra have been used to derive 
oxygen abundances from the triplet together with carbon abun- 
dances from the Ci lines in the HARPS spectra. Sect. 2 con- 
tains a description of the stellar spectra and the measurement 
of equivalent widths. The model-atmosphere analysis, includ- 
ing non-LTE corrections, is dealt with in Sect. 3, and the abun- 
dances obtained when adopting spectroscopic and photometric 
estimates of Teff and log g, respectively, are presented in Sect. 4 
and 5. Finally, Sect. 6 provides a comparison with other works 
and a discussion of the Galactic dispersion and evolution of C/O. 

2. Stellar spectra and equivalent widths 

Tables 1 and 2 in Delgado Mena et al. ( 120101 ) contain a list of 
100 dwarf stars with detected planets for which they have deter- 
mined C and O abundances based on spectra from the HARPS 
GTO survey. Out of this sample, 38 stars have FEROS spectra 
available in the ESO Science Archive. Of these, five stars have 
Teff around 5100K. The remaining 33 stars (listed in Table [TJ 
have 5400 < Teff < 6400 K, which is a suitable range for a dif- 
ferential abundance analysis with respect to the Sun. 

Pipeline extracted and wavelength calibrated spectra ob- 
served under the programmes listed in Table [T] were acquired 
from the ESO Science Archive. The HARPS spectra cover a 
wavelength range from 3800 to 6900 A and have a resolution of 
7? ^ 1 15 000. After combination of many individual spectra, the 
signal-to-noise {S IN) ratio ranges from 250 to more than 1000. 
The FEROS spectra range from 3500 to 9200 A with a resolution 
of :^48 000 . Most of these spectra have 5 /A^ from 200 to 300 
(in the region of the Oi triplet), except for HD 216770 that has 
SIN ^ 120. 

The spectra were normalized with the IRAF continuum task 
using a cubic spline fitting function of relatively low order, i.e. 
a wavelength scale > 20 A. The measurement of an equivalent 
width (EW) was carried out with reference to local continuum 



regions selected to be free of lines in the solar spectrum and 
situated within 4 A from the line. Care was taken to use the same 
continuum windows in all stars. 

The measurement of the equivalent widths of the O i triplet 
lines causes particular problems. These lines are somewhat 
strong and have very broad wings. If measured by fitting a Voigt 
profile, the equivalent width becomes very sensitive to the set- 
ting of the continuum. This is even a problem in solar spectra 
having extremely high resolution and S IN. Thus, the EWs of 
the O I triplet lines measured by Asplund et al. (2004 ) and Caffau 
et al. (120081) differ by 5 to 10 mA, which is a major reason for 
the difference in the solar oxygen abundance derived in these 
two 3D studies. In the present work, equivalent widths are, in- 
stead, measured by Gaussian fitting using the IRAF splot task, 
which for the resolution of FEROS turns out to provide nearly 
the same EW as obtained by direct integration over a range of 
±0.4 A centred on the line. The method neglects the contribution 
from the broad wings, but this was compensated for by also cal- 
culating the EW in the model-atmosphere analysis over a win- 
dow of ±0.4 A. The same procedure was followed for the two C i 
lines, although they are sufficiently weak to ensure that Gaussian 
fitting provides a good estimate of the total EW of the line. 

The measured equivalent widths are given in Table [1] It is 
noted that weak CN and C2 Unes are present in the region of 
the O I triplet, but fitting with a Gaussian profile minimizes the 
influence of such blends. In general, the equivalent widths of the 
O I lines have a precision of 2 - 3 mA and those of the C i lines a 
precision of about 1 mA. 

To make a reliable differential analysis with respect to the 
Sun, it is important to apply a solar flux spectrum observed in 
the same way as the stellar spectra. Fortunately, this is possi- 
ble; spectra of reflected sunlight from Ceres and Ganymede have 
been obtained with both HARPS and FEROS, and with S IN ra- 
tios similar to those of the programme stars. As seen from Table 
[T] the equivalent widths measured for the two minor planets 
agree very well. Therefore, the average values have been applied 
in the model-atmosphere analysis. 

3. Analysis 

3.1. Model atmospheres and line broadening 

For each star, a plane parallel (ID) model atmosphere was ob- 
tained from the standard MARCS grid (Gustafsson et al. 120081 ) 
by interpolating to the Teff , log g, and [Fe/H] values of the star, 
and the Uppsala program EQWIDTH was used to calculate 
equivalent widths of the C i and O i lines as a function of C or O 
abundance assuming LTE. Interpolation to the observed equiva- 
lent widths then yields the LTE abundances. 

Line data used in the analysis are given in Table |2] The gf 
values of the C i lines are based on quantum mechanical calcu- 
lations by Hibbert et al. (119931 ) and those of oxygen are taken 
from Hibbert et al. (I199H . Doppler broadening by non-thermal, 
small-scale motions in the atmosphere is described by one pa- 
rameter, ^tmb' in the usual way. For the Oi lines, collisional 
broadening caused by neutral hydrogen and helium atoms is 
based on the quantum mechanical calculations of Barklem et al. 
( 120001) . Their work does not include the two C i lines, for which 
the Unsold (1955) approximation with an enhancement factor of 
two was adopted as in Asplund et al. (i200 5 ). The Ci lines are 
sufficiently weak to be insensitive to the assumed enhancement 
factor; if it is changed from two to one, the derived [C/H] values 
change by less than 0.01 dex. 
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Table 2. Line data and deiived solar abundances. 



ID 


Wavelength 
(A) 


^exc 

(eV) 




(mA) 


A(X)o'' 
LTE 


A(X)o 
non-LTE 


Ci 


5052.17 


7.685 


-1.301 


35.9 


8.44 


8.43 


Ci 


5380.34 


7.685 


-1.616 


21.5 


8.44 


8.43 


Oi 


7771.94 


9.146 


0.369 


71.8 


8.88 


8.66 


Oi 


7774.17 


9.146 


0.223 


61.8 


8.86 


8.66 


Oi 


7775.39 


9.146 


0.002 


48.6 


8.84 


8.67 



Notes. For an element X, A(X) = log (A'x/A'h) + 12.0 



3.2. Non-LTE corrections 

It is well know that the strength of the Oi triplet is subject to 
large non-LTE effects (e.g. Kiselman l2001l) . The most advanced 
study is due to Fabbian et al. ('20091), who calculated non-LTE 
corrections for a grid of MARCS models ranging from 4500 to 
6500K in T^ff, from 2.0 to 5.0 in log g, and from -3.0 to 0.0 in 
[Fe/H]. Their model atom contains 54 energy levels and includes 
electron collision cross sections based on quantum mechanical 
calculations by Barklem ( I2007I I. An IDL program, made avail- 
able, allows interpolation of the non-LTE correction to the stellar 
values of T^ff, log g, and [Fe/H] for a given LTE oxygen abun- 
dance. For stars with [Fe/H] > 0.0, one needs to make an extrap- 
olation in [Fe/H], but because the correction changes only a little 
from [Fe/H] = -0.5 to 0.0 (see Fig. 7 in Fabbian et al.), this is 
not expected to introduce any significant errors. 

The "dark horse" in the non-LTE calculations is the han- 
dling of inelastic collisions with hydrogen atoms, because no 
experimental or quantum mechanical estimates of cross sections 
are available. Therefore, Fabbian et al. adopted the classical 
Drawin formula (Drawin,1968 ) scaled by an empirical factor S h- 
According to the 3D, non-LTE study of the solar centre-to-limb 
variation of the O i triplet lines by Pereira et al. (2009^ who ap- 
plied the same atomic model as Fabbian et al., = 0.85 leads 
to the best agreement with the observations. Because Fabbian 
et al. ( 120091 ) provide non-LTE coiTections for both Su - (H 
collisions are neglected) and 5h = 1 (Drawin's formula), it is 
possible to interpolate their corrections to 5h = 0.85, which is 
the value applied in the present paper 

According to Table |2] the solar non-LTE corrections range 
from -0.22 to -0.17 for the three O i lines, i.e. close to -0.20 dex 
on average. As seen from Figs. 7 and 8 in Fabbian et al. (2009), 
the coiTections become increasingly more negative for higher 
Jeff and lower \ogg, whereas stars with lower T^ff than the 
Sun have less negative non-LTE corrections. Altogether, this 
leads to significant differential non-LTE corrections of [O/H] 
ranging from about -0.2 dex at T^g ^ 6400 K to -i-0. 1 dex at 
Teff ^ 5400 K. 

The non-LTE corrections for the C i /15052 and /15380 lines 
are much smaller than those for the O i triplet. According to the 
detailed calculations of Takeda & Honda (2005) adopted in this 
paper, the corrections are -0.01 dex for the Sun and changes 
only slightly as a function of T^ff and log g. Their calculations 
were made on the assumption that S'h = L Assuming instead 
5h = 0, Asplund et al. (2005 ) find a solar non-LTE coiTection of 
-0.03 dex. In this case the coiTections of [C/H] would be some- 
what higher, but still close to be negligible compared to other 
error sources. 
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Fig. 1. Carbon abundances of Delgado Mena et al. (1201 Oi l in 
comparison with carbon abundances derived in this paper. 
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Fig. 2. Oxygen abundances of Delgado Mena et al. (120101) in 
comparison with non-LTE oxygen abundances derived in this 
paper. 



3.3. Random and systematical errors 

Random errors arise from the uncertainties of the measured 
equivalent widths and from the uncertainty of the adopted stel- 
lar atmospheric parameters. The first source is best estimated by 
comparing abundances obtained from different spectral lines of 
the same element. The [C/H] abundances derived from the two 
C I lines agree very well; the average standard deviation of the 
mean is cr[C/H] = 0.009 dex. From the three O i lines one gets 
cr[0/H] (non-LTE) = 0.017 dejfl Adding these errors in quadra- 
ture, the one-sigma error of [C/O] due to EW uncertainties be- 
comes ^ 0.02 dex. 

Based on the rms deviations between spectroscopic and pho- 
tometric parameters discussed in Sect. |5] the one-sigma errors 
for each set of parameters are estimated to be cr(reff) - 45 K, 
cr(logg) = 0.08, cr[Fe/H] = 0.04, and cr(^turb) = 0.07 km s"'. 
The effects on the derived [C/H], [O/H], and [C/O] values of 
the corresponding changes in stellar parameters are shown in 
Table|3] As seen, the largest contribution to the abundance errors 
arises from the error in Teff with [O/H] being more affected than 
[C/H]. This is due to a higher excitation potential of the O i lines 
(9.15eV) than in the case of the Ci lines (7.69 eV), and to the 
high Teff-sensitivity of the O i non-LTE corrections. Still, there 
is some canceling effect of the T^s dependence for the C i and 



- Interestingly, this number increases to 0.029 dex for the LTE values 
of [O/H] showing the importance of individual non-LTE corrections for 
the three O i lines. 
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Table 3. Changes in the derived non-LTE values" of [C/H], 
[O/H], and [C/O] due to the listed changes in Jeff, logg, [Fe/H], 
and^turb- 



Source 


A[C/H] 


A[0/H] 


A[C/0] 


Areii = +45 K 


-0.030 


-0.059 


+0.029 


Aloggr = +0.08 


+0.026 


+0.034 


-0.008 


A[Fe/H] = +0.04 


+0.000 


+0.001 


-0.001 


A^turb = +0.07 km S-' 


-0.002 


-0.005 


+0.003 



Notes. Calculated for a representative star, HD 134987 



the Oi lines; the error in [C/O] is only half the eiTor in [O/H]. 
Furthermore, the derived [C/H] and [O/H] values have nearly 
the same gravity dependence, so that the log gr-induced error on 
[C/O] becomes relatively small. Finally, it is noted that errors 
induced by the uncertainties in [Fe/H] and ^turb are negligible. 

When adding all eiTors in quadrature, the total eiTors become 
cr[C/H] = 0.04, cr[0/H] = 0.07, and cr[C/0] = 0.04. These 
are the random one-sigma errors, but systematical errors may 
also be present. As discussed in Sect. |5] there may be a system- 
atic error of Teff as a function of [Fe/H] ranging from +60 K at 
[Fe/H] = -0.5 to -60 K at [Fe/H] = +0.3. The coiTespond- 
ing change in [C/O] ranges from +0.04 at [Fe/H] = -0.5 to 
-0.04 at [Fe/H] - +0.3. There may also be an effect of apply- 
ing ID model atmospheres instead of more realistic 3D models. 
According to Asplund et al. (2004 ), the difference of the oxygen 
abundance derived from the O i triplet using a 3D and a MARCS 
model for the solar atmosphere is -0.06 dex, whereas the corre- 
sponding difference for the solar carbon abundance determined 
from the C i lines is -0.01 dex (Asplund et al. l2005l l. It remains 
to be seen if these 3D - MARCS differences change with T^ff or 
[Fe/H], but because the range of parameters for the present sam- 
ple of stars is less than ±600 K in Teff and +0.5 dex in [Fe/H], I 
do not expect any major 3D effects on the derived [C/O] val- 
ues. It is also noted that the uncertainty of the Su parameter 
for hydrogen collisions, say a change from the adopted value 
5h = 0.85 to an extreme of =0, has only a small effect on 
[C/O]; the changes are within ±0.02 dex. 

4. Results based on spectroscopic stellar 
parameters 

Carbon and oxygen abundances are first derived with MARCS 
model atmospheres having the same parameters as applied by 
Delgado Mena et al. ( 120 1 Ob . The original source of these param- 
eters is Sousa et al. (2008), who in a differential, LTE analysis 
with respect to the Sun determined Teff, log g, [Fe/H], and ^turb 
by requesting that iron abundances derived from Fe i and Fe ii 
lines have no systematic dependence on excitation potential, ion- 
ization stage, or equivalent width. The parameters are listed in 
Table ID together with the LTE and non-LTE values of [C/H] and 
[O/H] derived in this paper 

As seen from Fig. [1] there is an excellent agreement be- 
tween the carbon abundances derived by Delgado Mena et al. 
(1201 Oi l and the values derived in this paper The rms deviation in 
[C/H] is only 0.025 dex. This is not surprising, because in both 
studies C abundances are derived from equivalent widths of the 
.i.i5052, 5380 C I fines in fiigfi S /N HARPS spectra using tfie 
same set of stellar parameters. 

Tfie situation is, fiowever, very different wfien one compares 
tfie oxygen abundances derived by Delgado Mena et al. (,2010) 
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from tfie /16300 [O i] line witfi O abundances derived in tfiis pa- 
per from FEROS spectra of tfie AlllA O i triplet (see Fig. |2]l. In 
tfiis case, tfie rms deviation is as large as 0.12 dex. Tfie prob- 
lem can be illustrated by comparing tfie spectra of two stars 
witfi similar atmospfieric parameters Teff, log g, and [Fe/H], i.e 
HD 196050 (5917 K, 4.32, 0.23) and HD 216435 (6008 K, 4.20, 
0.24). As seen from Figs. |3]and|4] tfie Ci and Oi lines fiave 
about tfie same strengtfis in tfie two stars. Actually, tfie lines 
of HD 216435 are a bit wider and fiave sligfitly larger EWs, 
tfian tfiose of HD 196050, but due to tfie small differences in 
Teff and log g between tfie two stars, tfiis is to be expected for 
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fixed values of [C/H] and [O/H]. As seen from Table |4] the 
model-atmosphere analysis provide very similar O abundances, 
i.e. [O/H] = 0.12 for HD 196050 and [O/H] = 0.09 in the case 
of HD 216435. This is in stark contrast to the results of Delgado 
Mena et al. ( l20T0b . who obtain [O/H] = 0.36 for HD 196050 and 
[O/H] = -0.02 for HD 216435. 

Figure |5] shows a comparison of HARPS spectra of 
HD 1 96050 and HD 2 1 6435 in the 6300 A region. The [O i] +Ni i 
blend is affected by a telluric O2 line in many of the available 
HARPS spectra of HD 196050, so the spectrum shown in Fig. 
|5]is an average of a subset of spectra from July 2004, when the 
topocentric radial velocity of HD 196050 shifts the telluric line 
away from the [O i] line. As seen from the figure the equivalent 
width of the [Oi]+Nii blend is nearly the same in the spectra 
of HD 196050 and HD 216435, i.e. EW = 8.0 and 7.5 mA, re- 
spectively. Adopting [Ni/H] = 0.28 and 0.29 for the two stars 
(as used by Delgado Mena et al.) and solar flux data ^^[Oi] - 
3.5 mA and ^^(Nii) = 2.0mA (Caffau et al. 2008 ), a differen- 
tial analysis with respect to the Sun leads to [O/H] - 0.17 + 0.10 
for HD 196050 and [O/H] = 0.1 1 + 0.10 for HD 216435, where 
the main contribution to the error arises from the uncertainty of 
setting the continuum level. These [O/H] values agree well with 
the non-LTE values derived from the O i triplet, but lend no sup- 
port to the large difference in [O/H] found by Delgado Mena 
et al. (120101 ). Perhaps, they applied a different set of HARPS 
spectra for HD 196050; in any case, the analysis shows that oxy- 
gen abundances derived from the [Oi]+Nii blend in spectra of 
metal-rich stars are uncertain due to difficulties in setting the 
continuum level and correcting for the large contribution of the 
Ni I line. 

5. Results based on photometric stellar parameters 

In order to test the reliability of the derived C/O ratios, pho- 
tometric atmospheric parameters of the stars have been deter- 
mined as an alternative to the spectroscopic parameters adopted 
by Delgado Mena et al. (120101) . Teff was taken as the mean value 
derived from the (b-y) and (V-K) colour indices using the 
IRFM calibrations of Casagrande et al. (2010) with V magni- 
tudes and (b-y) from Olsen ( I1983I I and K magnitudes from the 
2MASS catalogue (Skrutskie et al. I2006i All stars are closer 
than 60 pc according to their Hipparcos parallaxes, so reddening 
corrections can be neglected. 

The surface gravities of the stars were determined from the 
fundamental relation 

a M T rr 

log ^ = log — + 4 log -A + 0.4(Mboi - MboLo) (1) 

where M is the mass of the star and Mboi the absolute bolo- 
metric magnitude. Hipparcos parallaxes (van Leeuwen 12007) 
were used to derive My and the bolometric correction adopted 
from Casagrande et al. (1201 Oi l. The stellar mass was obtained 
by interpolating in the luminosity - log Jeff diagram between the 
Yonsei -Yale evolutionary tracks by Yi et al. (i2003j as described 
in detail by Nissen & Schuster (2012]l. Due to the small error 
of the Hipparcos parallaxes, the estimated error of log g is only 
0.04 dex. 

The metallicity was determined from 12 Feii lines with 
equivalent widths being measured from the HARPS spectra. By 
using Fe 11 lines only, the derived [Fe/H] value is unlikely to be 
affected by departures from LTE (Mashonkina et al. 1201 11 Lind 
et al. 120121) . The Feii lines have equivalent widths in the range 
10 i EW i 90 mA, which allows a determination of ^tuit from 
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Fig. 6. Differences between photometric gravities determined 
via Hipparcos parallaxes and spectroscopic gravities used by 
Delgado Mena et al. (l20T0t 
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Fig. 7. Differences between photometric values of T^s deter- 
mined from colour indices and spectroscopic values used by 
Delgado Mena et al. (12010b . 
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Fig. 8. Difference between [C/O] values derived in this paper, 
when using photometric and spectroscopic stellar parameters, 
respectively. 

the requirement that the derived [Fe/H] should not depend on 
EW. 

It is noted that this procedure for the determination of stellar 
parameters has to be iterated, because the T^ff calibrations, bolo- 
metric corrections and mass determinations depend on [Fe/H]. 

When comparing the photometric and spectroscopic param- 
eters, the following rms deviations are obtained: rms(Areff) = 
62 K, rms(Alog g) = 0.12, rms(A[Fe/H]) = 0.05, and rms(A^tui-b) 
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= O.lOkms"'. The scatter for [Fe/H] and ^turb is as small as one 
would expect, if the one-sigma errors of each of the two studies 
are cr[Fe/H] ^ 0.03 and o-{^t^rh) - 0.07 km s~^ In the case of 
log g, the scatter is larger than expected from the estimated er- 
rors, i.e. cr(log g) - 0.06 for the spectroscopic gravities (Sousa 
et al. 120081) and cr(log g) = 0.04 for the photometric values ob- 
tained via the Hipparcos parallaxes. As shown in Fig. |6] this is 
mainly due to a significant dependence of Alog g on Teff ; the rms 
scatter around the fitted line is only 0.07 dex. This trend is not 
easy to explain; according to Lind et al. ( 120121 Fig. 4) non-LTE 
effects on the ionization balance of Fe (used to determine the 
spectroscopic gravities) are almost negligible for the effective 
temperatures, gravities and metallicities of the present sample of 
stars. 

A similar problem is present, when one compares photomet- 
ric and spectroscopic temperatures. In this case, the difference 
depends on [Fe/H] as shown in Fig. [7] Here, the rms deviation 
around the fit is 40 K. Again it is unclear, what is causing this 
trend with [Fe/H]. 

Despite of these systematic differences between the photo- 
metric and spectroscopic stellar parameters, the C/0 ratios de- 
rived in this paper are not much affected as seen from Fig. |8] 
[C/O] is immune to changes in log g, whereas the trend of AT^eif 
with [Fe/H] shown in Fig. |7] introduces a small dependence of 
A[C/0] on [Fe/H]. 

6. Discussion 

Fig.|9]shows the C/O ratios derived in this paper compared to the 
ratios derived by Delgado Mena et al. (i2010j with the same set 
of spectroscopic stellar parameters applied in the two studies. 
Evidently, the C/O ratios of Delgado Mena et al. have a much 
larger scatter than the ratios obtained in this paper Furthermore, 
the Delgado Mena et al. C/O ratios tend to be higher than the 
values in the present study, which is mainly due to the fact that 
they adopt a solar ratio (C/O)© = 0.66 instead of the value 
(C/O)0 = 0.58 corresponding to the solar abundances derived 
in this paper (see Table HJ. Recent 3D studies lead to nearly the 
sa me C/O ratio in the Sun, i . e (C/0 )q - 0.55 + 0.10 (Asplund et 
al. 120091 Caffau et aL lIOMHIUTU I). The solar C/O ratio adopted 
by Delgado Mena et al. is on the high side of the quoted error 
range for these 3D studies. 

Excluding two stars with [Fe/H] < -0.3 having thick-disk 
kinematics, one obtains a Hnearfit, C/O = 0.58 -i- 0.48 [Fe/H], to 
the data of this pape. If the two thick-disk are included the rela- 
tion becomes C/O = 0.56 + 0.54 [Fe/H]. In both cases, the rms 
dispersion around the fit is cr (C/O) = 0.06. In comparison, the 
Delgado Mena et al. (120 1 01) data have a scatter of cr = 0.14 rel- 
ative to a linear fit with three stars having C/O - 1 and 10 stars 
having C/O > 0.8. As discussed in Sect. H) this difference in the 
C/O distribution is mainly due to the large differences in the oxy- 
gen abundances derived from the /J6300 [O i] line by Delgado 
Mena et al. and the O abundances derived in this paper from 
the Oi triplet as exemplified by the comparison of HD 196050 
and HD 215435. The small dispersion of C/O derived in this pa- 
per suggests that the high C/O values (> 0.8) found by Delgado 
Mena et al. (2010) and by Petigura & Marcy (2011) are spuri- 
ous due to difficulties in determining oxygen abundances from 
the /16300 [O i]-i-Ni i blend especially at high metallicities. In the 
solar spectrum, the EWs of the Ni i and [O i] lines are 2.0 and 
3.5 mA, respectively, (Caffau et al. 2008) i.e. a ratio of R ^ 0.55, 
but in a solar-type star with [Fe/H] - +0.3, [O/Fe] = -0.15, and 
[Ni/Fe] = +0. 1 (Bensby et al. 120051 . the ratio is increased to 
^ ^ 1. 
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Fig. 9. Carbon-to-oxygen ratios derived in this paper (filled, 
black circles) in comparison with the C/O ratios determined by 
Delgado Mena et al. (2010) (open, red circles). The error bars 
shown refer to data from this paper In the case of HD 196050 
and HD 216435, the filled and open circles are connected by a 
straight line. 



In view of the possible systematic errors in the 3D, non- 
LTE corrections for the C and O abundances derived from high- 
excitation C I and O i lines, it is important to check the trend and 
scatter of C/O obtained in this paper by independent determina- 
tions of C/O. Such a study has been carried out by Bensby et al. 
( 120041 and Bensby & Feltzing ( f2006l ). who determined O and C 
abundances from weak forbidden lines, i.e. [Oi] at 6300 A and 
[C i] at 8727 A. The non-LTE corrections for these lines are van- 
ishing and 3D corrections are small in the solar case (Asplund et 
al. 2004, 2005; Caffau et al. 2008, 2010), but the Hnes are weak 
and affected by blends; in addition to the Ni i blend of the /16300 
[Oi] line, a weak Fei line is overlapping the /18727 [Ci] line. 
Hence, very high resolution and high S /N spectra are needed 
like those obtained with the Coude Echelle Spectrograph at the 
ESO 3.6 m telescope by Bensby et al. 

Figure [10] shows [C/O] versus [Fe/H] for a sample of 42 
dwarf stars with 5000 < T^ff < 6500 K from Table 1 in Bensby 
& Feltzing (2006), divided into thin- and thick-disk stars ac- 
cording to their kinematics. In addition, the figure shows data 
from the present study based on model atmospheres with pho- 
tometric parameters. For this sample, three stars with a total 
space velocity relative to the Local Standard of Rest Vlsr > 
80kms '(Holmberg et al. 120091 ) have been classified as belong- 
ing to the thick disk. Four stars are in common between the two 
samples; they show good agreement in [C/O], i.e. a mean and 
rms deviation (Bensby - this paper) A[C/0] = 0.02 ± 0.05. 

As seen from Fig. [10] there is an excellent agreement be- 
tween [C/O] data from Bensby & Feltzing (2006 ) and from this 
paper, including a systematic difference in [C/O] between thin- 
and thick-disk stars at low metallicity. A linear fit to the thin-disk 
data yields a trend 



[C/O] = -0.002 + 0.22 [Fe/H], 



(2) 



with a standard deviation of +0.048 dex for 30 stars from this 
work and +0.068 dex for 27 stars from Bensby & Feltzing. This 
equation corresponds to [C/O] - 0.086 at [Fe/H] = +0.4 or C/O 
= 0.71, when adopting a solar ratio (C/O)© - 0.58. If the solar 
ratio is (C/O)© = 0.65, Eq. (2) corresponds to C/O = 0.79 at 
[Fe/H] = +0.4. 
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Fig. 10. [C/O] versus [Fe/H]. Data from this paper are plotted 
with circles and data from Bensby & Feltzing (2006') are shown 
with (red) triangles. Open symbols refer to stars with thin-disk 
kinematics and filled symbols to thick-disk stars. A linear fit to 
the thin-disk stars, given in Eq. (2), is shown by a dashed line. 
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Fig. 11. [C/O] versus [Fe/H] with data from Takeda & Honda 
(120051 1 (circles) and Ramirez et al. (|2009) (red triangles). Eq. (2) 
is shown as a dashed line. 



Another two high-precision studies of carbon and oxygen 
abundances in solar-type stars (Takeda & Honda 12005 1 Ramirez 
et al. I2009I I also show good agreement with the [C/0]-[Fe/H] 
trend found in this paper as seen from Fig. [TT] In both works, the 
adopted abundances are obtained from high-excitation C i and 
O I lines with non-LTE corrections based on their own calcula- 
tions. The [C/O] values from Takeda & Honda, which include 
27 planet-harbouring stars, distribute around the line given in 
Eq. (2) with an rms deviation of +0.070 dex, and the data from 
Ramirez et al. for solar twins and analogs have an rms deviation 
of 0.058 dex. In both cases, the scatter can be explained by the 
errors of the abundance determination, and although a few stars 
from Takeda & Honda fall above [C/O] - 0.14, which corre- 
sponds to C/O - 0.8, this may be accidental. 

While oxygen is probably formed on a relatively short time- 
scale in massive stars exploding as Type II SNe, the origin of 
carbon is more uncertain with many possible sources including 
Type II supernovae, novae, winds from massive stars, and low- 
to intermediate-mass stars (Gustafsson et al. 11999) . As shown 



by Akerman et al. (2004), [C/O] raises from a level of -0.5 dex 
in halo stars with [O/H] < -1 to [C/O] ^ in disk stars with 
solar metallicity. A similar trend is found from H ii regions in 
spiral and irregular Galaxies (Garnett et al. 120041 . From this ap- 
parent universal metallicity dependence, Akerman et al. argue 
that the carbon enrichment in the Galactic disk is mainly due 
to metal-dependent winds from massive stars. In order to ex- 
plain the difference in [C/O] between thin- and thick-disk stars 
shown in Fig. (TO] it is, however, necessary to invoke a signifi- 
cant carbon contribution from low- and intermediate-mass stars 
ocurring on a relatively long time-scale in the Galactic thin disk 
(Chiappini et al. 2003. Carigi et al. l2005l Cescutti et al. ,2009l l, 
whereas data for [C/O] in stars belonging to the thick disk and 
the Galactic bulge may be explained by metallicity-dependent C 
and O yields for massive stars alone (Cescutti et al. 2009i. The 
situation is, however, complicated and the [C/O] data at high 
metallicity are not well explained, so more precise observations 
of C/O, more reliable yields, and realistic Galactic models are 
needed to advance the field. 

7. Conclusions 

The results in this paper and those of Bensby & Feltzing (I2006I I. 
Takeda & Honda (2005), and Ramirez et al. ( 120091 ) show a tight 
[C/0]-[Fe/H] correlation for thin-disk stars. Furthermore, [C/O] 
is not reaching higher than [C/O] ^0.1 dex (corresponding to 
C/O ^ 0.7 for a solar ratio of (C/0)o = 0.55) even at the highest 
metallicities. This suggests that the recent findings by Delgado 
Mena et al. (2010) and Petigura & Marcy (201 1 ) of C/O > 0.8 in 
a significant fraction (10 - 15 %) of solar-type stars are spurious 
due to difficulties in measuring the strength of the weak /16300 
[Oi]+Nii blend with sufficient precision and to correct for the 
Ni I contribution when deriving oxygen abundances. 

As mentioned in the introduction, modelling of planet for- 
mation (Bond et al. 1201 Oi l indicates ffiat C/O > 0.8 in proto- 
planetary disks with uniform composition lead to the formation 
of terrestrial "carbon planets" consisting of carbides and graphite 
instead of "Earth-like" silicates. Assuming that the composition 
of a proto-planetary disk is the same as that of the host star, the 
results in this paper (C/O < 0.8) lend no support to the existence 
of such carbon planets. 

The separation of [C/O] between thin- and thick-disk stars 
suggests that carbon in the thin disk are mainly formed in low- 
and intermediate-mass stars. Assuming that oxygen comes from 
massive stars exploding as Type II SNe, the small dispersion of 
[C/O] indicates that the nucleosynthesis products of these ob- 
jects are well mixed in the interstellar medium before new stars 
are formed. A detailed understanding of the formation of the two 
elements requires, however, further observational and theoretical 
studies. 
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Table 1. ESO observing programme numbers, SjN ratios, and equivalent widths measured by Gaussian fitting to line profiles. 
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Oi Mil 5 


ID 


Programme 


SIN 


Programme 




EWimk) 


EW{mk) 


EW{mk) 


EWimk) 


EW{mk) 


HD 142 


72.C-0488 


450 


83.A-9011 


300 


69.8 


45.2 


131.7 


114.3 


95.8 


HD1237 




350 


60.A-9700 


200 


30.9 


16.8 


59.1 


50.2 


36.4 


HD4308 




600 


74.D-0086 


500 


25.5 


14.1 


68.1 


58.5 


42.4 


HD 16141 




300 


83.A-9011 


250 


45.2 


27.9 


86.3 


75.7 


58.6 


HD 20782 




900 




200 


33.5 


19.9 


72.2 


59.5 


49.3 


HD 23079 
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84.A-9004 


200 


34.1 


19.7 


80.8 


68.5 


51.8 
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250 
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46.2 
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250 
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HD 108147 
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55.6 


35.0 


114.9 


98.1 


80.9 


HD 111232 




550 
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19.6 


10.5 


54.8 


45.5 


33.0 


HD 1 14729 




900 
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31.6 


18.2 


76.0 


64.9 


52.6 


HD 117618 




600 




300 


44.8 


28.1 


92.1 


79.4 


60.2 


HD 134987 
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33.8 
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950 


83.A-9011 


300 


56.8 


37.1 


97.5 


82.6 


69.5 


HD 202206 




700 




200 


43.3 


26.7 


74.7 


67.0 


48.7 


HD 210277 




900 




200 


41.1 


25.0 


65.1 


57.7 


44.4 


HD 212301 


82.C-0312 


350 


85.C-0743 


350 


58.4 


37.2 


112.8 


97.6 


81.8 


HD 213240 


72.C-0488 


300 


83.A-9011 


200 


51.9 


33.5 


99.1 


82.5 


71.1 


HD 216435 




700 




200 


61.2 


40.1 


108.9 


96.5 


77.1 


HD 216437 


80.D-0408 


700 




200 


36.7 


21.6 


52.4 


47.3 


37.3 


HD 216770 


72.C-0488 


600 




120 


54.9 


36.9 


96.3 


81.8 


65.7 


Ceres 


60.A-9036 


400 


85.A-9027 


350 


36.3 


21.4 


72.3 


61.2 


48.2 


Ganymede 




450 


77.C-0766 


250 


35.6 


21.6 


71.3 


62.4 


49.0 
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Table 4. Adopted atmospheric parameters and derived carbon and oxygen abundances 



ID 


^eff 


logs' 




j- 

Sturb 


LC/HJ 


LC/HJ 


LO/HJ 


LO/HJ 


c/u 




(K) 






kms~* 


LTE 


non-LTE 


LTE 


non-LTE 


non-LTE 


HD142 


6403 


4.62 


0.09 


1.74 


0.21 


0.21 


0.21 


0.08 


0.78 


HD 1237 


5514 


4.50 


0.07 


1.09 


0.07 


0.07 


0.07 


0.14 


0.50 


HD4308 


5644 


4.38 


-0.34 


0.90 


-0.15 


-0.15 


0.04 


0.09 


0.34 


HD 16141 


5806 


4.19 


0.16 


1.11 


0.07 


0.06 


0.12 


0.06 


0.59 


HD 20782 


5774 


4.37 


-0.06 


1.00 


-0.06 


-0.07 


-0.03 


-0.03 


0.54 


HD 23079 


5980 


4.48 


-0.12 


1.12 


-0.14 


-0.14 


-0.09 


-0.12 


0.55 


HD28185 


5667 


4.42 


0.21 


0.94 


0.22 


0.22 


0.13 


0.14 


0.71 


HD 30177 


5588 


4.29 


0.38 


1.08 


0.38 


0.38 


0.31 


0.31 


0.68 


HD 52265 


6136 


4.36 


0.21 


1.32 


0.15 


0.15 


0.18 


0.05 


0.72 


HD 65216 


5612 


4.44 


-0.17 


0.78 


-0.20 


-0.19 


-0.09 


-0.04 


0.40 


HD 69830 


5402 


4.40 


-0.06 


0.80 


-0.05 


-0.05 


-0.10 


-0.02 


0.55 


HD 73256 


5526 


4.42 


0.23 


1.11 


0.22 


0.22 


0.13 


0.18 


0.64 


HD 75289 


6161 


4.37 


0.30 


1.29 


0.11 


0.10 


0.14 


0.00 


0.73 


HD 82943 


5989 


4.43 


0.26 


1.10 


0.20 


0.20 


0.22 


0.14 


0.66 


HD 92788 


5744 


4.39 


0.27 


0.95 


0.24 


0.24 


0.21 


0.18 


0.66 


HD 108147 


6260 


4.47 


0.18 


1.30 


0.06 


0.06 


0.13 


0.00 


0.66 


HD 111232 


5460 


4.43 


-0.43 


0.62 


-0.17 


-0.16 


0.04 


0.13 


0.30 


HD 114729 


5844 


4.19 


-0.28 


1.23 


-0.21 


-0.21 


-0.09 


-0.12 


0.47 


HD 117618 


5990 


4.41 


0.03 


1.13 


0.03 


0.03 


0.05 


-0.01 


0.63 


HD 134987 


5740 


4.30 


0.25 


1.08 


0.29 


0.29 


0.24 


0.20 


0.70 


HD 160691 


5780 


4.27 


0.30 


1.09 


0.27 


0.27 


0.25 


0.19 


0.69 


HD 168443 


5617 


4.22 


0.06 


1.21 


0.14 


0.14 


0.15 


0.15 


0.56 


HD 169830 


6361 


4.21 


0.18 


1.56 


0.13 


0.12 


0.29 


0.07 


0.66 


HD 179949 


6287 


4.54 


0.21 


1.36 


0.13 


0.12 


0.16 


0.04 


0.70 


HD 183263 


5991 


4.38 


0.34 


1.23 


0.22 


0.22 


0.15 


0.06 


0.83 


HD 196050 


5917 


4.32 


0.23 


1.21 


0.23 


0.23 


0.20 


0.12 


0.75 


HD 202206 


5757 


4.47 


0.29 


1.01 


0.16 


0.16 


0.10 


0.08 


0.70 


HD 210277 


5505 


4.30 


0.18 


0.86 


0.25 


0.25 


0.20 


0.24 


0.59 


HD 212301 


6271 


4.55 


0.18 


1.29 


0.12 


0.12 


0.13 


0.02 


0.73 


HD 213240 


5982 


4.27 


0.14 


1.25 


0.10 


0.10 


0.12 


0.03 


0.68 


HD 216435 


6008 


4.20 


0.24 


1.34 


0.20 


0.19 


0.22 


0.09 


0.73 


HD 216437 


5887 


4.30 


0.25 


1.31 


0.23 


0.22 


0.18 


0.11 


0.75 


HD 2 16770 


5424 


4.38 


0.24 


0.91 


0.25 


0.25 


0.12 


0.18 


0.68 


Sun 


5777 


4.44 


0.00 


1 .00 


0.00 


0.00 


0.00 


0.00 


0.58 



